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Abstract 

In this work we study the cosmological evolution of a dark energy model with two scalar 
fields, i.e. the tachyon and the phantom tachyon. This model enables the equation of state 
w to change from w > — 1 to w < — 1 in the evolution of the universe. The phase-space 
analysis for such a system with inverse square potentials shows that there exists a unique 
stable critical point, which has power-law solutions. In this paper, we also study another 
form of tachyon-quintom model with two fields, which voluntarily involves the interactions 
between both fields. 
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I. INTRODUCTION 



Recent observational data [H, 0, [sj strongly indicate that the Universe is spatially flat 
and accelerating at the present time. Within the framework of general relativity, cosmic 
acceleration can be sourced by an energy-momentum tensor which has a large negative 
pressure called dark energy (See Ref. Q for a recent review). The simplest candidate for 
dark energy seems to be a small positive cosmo logical constant, but it suffers from difficulties 
associated with the fine tuning and coincidence problem. This problem can be alleviated 
in models of dynamically evolving dark energy called quintessence jsj], which have tracker 
like properties where the energy density in the fields track those of the background energy 
density before dominating today. The phantom, whose kinetic energy term has the reverse 
sign, has been also proposed as a candidate of dynamical dark energy [fL There has been 
the enormous variety of DE models suggested in the literature, see Ref. 7[ for reviews. 

The analysis of the properties of dark energy from recent observations mildly favor models 
with lo crossing —1 in the near past 0,0]. But, neither quintessence nor phantom can fulfill 
this transition. The quintom scenario of dark energy is designed to understand the nature of 
dark energy with u across —1. The first model of quintom scenario of dark energy is given 
in Ref. [91] with two scalar fields, where one is quintessence and the other is phantom. This 
model has been studied in detail later on 10, 11, 12, 0, Q, 0], and recently a new type 
of quintom model inspired by the string theory has also been proposed, which only have a 
single scalar field 
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The role of the rolling tachyon [17[ in string theory has been widely studied in cos- 
mology, see Refs.[18l. Il9| . and especially Refs.|2C|, |2l|, |22j, |23( for dark energy. Some sort 
of tachyon condensate may described by effective field theory with a Lagrangian density 
C = — V((f))^/l + g^d^cpducf). It can act as a source of dark matter or inflation field. Mean- 
while the tachyon can also act as a source of dark energy depending upon the form of the 
tachyon potential. However, as compared to canonical quintessence, tachyon models require 
more fine-tuning to agree with observations. In Ref. 24|, the authors consider the Born- 
Infeld type Lagrangian with negative kinetic energy term. The Lagrangian density they 
choose is C = —V{ip)yJ\ — g^ v d^ipd v ip. It is clear that for the spatially homogeneous scalar 
field, the equation of state to = — 1 — 2 will be less than —1 unless the kinetic energy term 



0. This field is called phantom tachyon, see for example 25 



In principle, we can consider the multi-fields model including multiple tachyon and mul- 
tiple phantom tachyon. However, without loss of generality, we only consider the case of one 
tachyon and one phantom tachyon, since it is the simplest one to realize that the equation of 
state w cross —1 during the evolution of the universe and it shows most of the central ideas 
of such model. Because of the quintom-like behavior it shows, we call it tachyon-quintom 
for the convenience. This paper is organized as follows: in section II we study in detail 
the tachyon-quintom model with inverse square potentials. The numerical analysis shows 
this model is not sensitive to the initial kinetic energy density of tachyon and phantom 
tachyon, and we give the reason in detail. Then the phase-space analysis of the spatially 
flat FRW models shows that there exist a unique stable critical point, and we compare it 
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with tachyon model; in section III we present another two-field model which include the 
interaction between two fields; the section IV is summary. 



II. THE TACHYON-QUINTOM MODEL 

A. The tachyon-quintom model 

We assume a four-dimensional, spatially-flat Friedmann-Robertson- Walker Universe filled 
by a homogeneous tachyon with potential V(<f>), a homogeneous phantom tachyon <p with 
potential V(tp) and a fluid with barotropic equation of state p 7 = (7 — l)p 7 , < 7 < 2, such 
as radiation (7 = 4/3)or dust matter (7 = 1). In this section, we turn our attention to the 
possibility of the tachyon and phantom tachyon as a source of the dark energy. 

The action for such a system is given by 

(1) 



r , i (Mpi \ 

= d xy/-g I — h £</, + C v + C m I 



where M p is the reduced Planck mass, 1Z is the scalar curvature, C m represents the 
Lagrangian density of matter fields and 

£ = -V{4>) y/l + g^d^cfr (2) 

£^-%)Vl-j"V # (3) 

We now restrict to spatially homogeneous time dependent solutions for which <9j0 = di<p = 
0. Thus the energy densities and the pressure of the field and (p are given, respectively, by 

v{<f>) 



P<t> = r—— . P^ = "U(0)Vl-0 2 , (4) 

V 1 - 2 

^ = 4==> P v = -V(<pWl + <P l (5) 
V 1 + f 

Here a dot is derivation with respect to synchronous time. The background equations of 
motion are 

'^ = (6) 



1-0 2 v(<j>) d<p 

* 3/7^- 1 ^M=0 (7) 



1+0 2 V(^) ^ 

p 7 = -^H Pl (8) 



^ = -2Mt(P<P + P<t> + P<t + P<t + Pi + Pi) 



2M„ 2 



g 2 v(y) ■ ^ 
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together with a constraint equation for the Hubble parameter: 



. 1 V(t) V{f) 

H to + 7TT? + ftl (10) 



The potentials we considered are inverse square potentials: 



V^) = Ml<\>-\ V^) = M^- 2 (11) 



Those potentials allow constructing a autonomous system [26j, |27| using the evolution 



equations, and give power-law solutions. The cosmological dynamics of the tachyon field 



with inverse square potential was studied in Ref. [20, 122, l22|. Interestingly, the inverse 
square potential plays the same role for tachyon fields as the exponential potential does for 
standard scalar fields. 

We define the following dimensionless quantities : 

x * V4> - 3H*Mf * Vip ~ 3# 2 Mf z ~ 3# 2 m 2 1 ' 

Now the Eqs. (jTUl) and (Q can be rewrite as follow: 

1 - [l " ' ' h +z = n DE + z (13) 



H' 3 / y^-xl) yvil + x 2 ) 

H = -.\-^=r - ^frw +l] (14) 

where Qde measure the dark energy density as a fraction of the critical density, a prime 
denotes a derivative with respect to the logarithm of the scale factor, N = In a. 
Then the evolution Eqs. (jSJ) and (jjTj) can be written to an autonomous system: 




.r. 



-3(a^- y/Pfy^l - xl) (15) 



= 3y<t> — / 7= : g - Vp<i>y<f> x <f> + 7 



(16) 



-3(a: v + V^)(l + 4) (17) 



, yja - 4) y y (7 + 4) /?5 — 1 M , M 
Zy v I ! 7= ; I - VP<py<P x v + 7 I (18) 



1 — X 2 l/ l tl 



where 



4M 2 _ 4M 2 

3M 2 ' ^~3Mj (19) 



The equation of state of the dark energy is 



x- 
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FIG. 1: Evolution of the equation of state (cj)and density parameters'^) as a function of 
N for the dark energy model with = /?„ = 1/3, 7 = 1. Initial conditions (at N = —8) 
: a. solid line: x<$n = 0.9999999, ^ = 6 x 1CT 11 , x^i = 2.0, y v i = 6.5 x 10~ n ; b. dashed 
line: x^ = 1 x 10~ 12 , y^ = 6 x 10~ n , x V i = lx 10 -12 , y v i = 6.5 x 10~ n ; c. dotted line: x^ = 
0.5, y^i = 6 X 10~ n , x^i = 1.0, y^i = 6.5 x 10~ n . 

B. Numerical analysis 

Mapping between the number of e-foldings and the redshift z = a /a—l = 1/a—l, we note 
that at big bang nucleosynthesis (BBN) N-qbn ~ — 20 (z ~ 10 9 ), at matter-radiation equality 
N eq ~ — 8 (2 ~ 3200). We choose iV = — 8 as the initial number of e- folds, so choose the 
7 = 1 in Eqs. ffTBT) and (fT8l) is a good approximation. The evolutions of u and Qde are shown 



in Fig.HJ In Ref. 28| . the authors use standard Big Bang Nucleosynthesis and the observed 
abundances of primordial nuclides to give a constraints on the scalar matter : Qde < 0.045, 
at temperatures near 1 MeV. The initial values of x^^y^^x^ and y v given below are safely 
satisfy this requirement, since Opg < 6 x 10~ 7 at iV = —8, and the energy density of and (p 
are decreasing more slowly than the fluid (u^ = y. = _i-|_0 2 < 0,^ = y- = —l — <p 2 < — 1). 
So at temperatures near 1 MeV, the QLde is smaller. 

From pictures we can see that this model is not sensitive to the initial kinetic energy 
density of tachyon and phantom tachyon — <j>, x v = ip). When x^ = 0.9999999, y^ = 
6 x 10~ n , x^ = 2.0, y vi = 6.5 x 10 -11 , the current u and Q DE are -1.025896 , 0.722306, 
respectively. When i^ = lx 10 -12 , = 6 x 10~ n , x vi = 1 x 10~ 12 , y vi = 6.5 x 10~ n , 
the current u and Qde are —1.025841 , 0.722367, respectively. From Eq. (jlj), we know that 
the initial energy density of the tachyon varied by nearly four orders of magnitude is still 
consistent with current observational constraints. But the initial potential energy density of 
tachyon and phantom tachyon require fine-tuning to agree with observations. 

At the present, we want to explain in rough detail how solutions converge to the common 
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FIG. 2: Evolution of the x^, ytj), Xip and y^ as a function of N for the dark energy model with 
P<j> = /3u> = 1/3, 7 = 1. Initial conditions (at N = —8) : a. solid line: x^i = 0.9999999, y^ = 
6 x 1CT 11 , x V i = 2.0, y v i = 6.5 x 10~ n ; b. dashed line: x^ = 1 x 10~ 12 , y^ = 6 x 10~ n , x V i = 
1 x 10~ 12 , y v i = 6.5 x 10 -11 ; c. dotted line: = 0.5, = 6 x 10~ n , x V i = 1.0, y v i = 6.5 x 10~ n . 

solution for the different initial conditions which given in Fig.[TJ (In next subsection, we will 
know that there is only one stable critical point, so the solutions converge to a common, 
cosmic evolutionary track is not surprised. ) 

From the Fig.0 we can see that the evolution of the x^, y^, x v and function 
of N. We observe that the and y v evolutionary tracks are not sensitive to the initial 
conditions of x^ and x v which given in Fig.[TJ This can be seen from the Eqs. (|T6|) and 
( TTSjl . The initial value of x^, y^, x v and y v are so small that we can safely neglect the 
first three term of the right hand side of Eqs. f|T6|) and ([TBI . So when —8 < N < —2, we 
get y'^ ~ 37?/^, y' v ~ 37?/^. When N > —2 , the different initial value of x^ and x v have 
converged to a common evolutionary track, so from then on even the value of x^, y^, x^ and 
y v become large, the convergence of y^ and y v with different initial value of x^ and x v is no 
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change. 

Now, consider the case (a) in which the initial value of x$ a is large, such as the solid 
line in Fig.[2] . Since x^ > yj fyy^ , so x^ a will decrease until x^ = ^ fyy^a = x A (see 
Eq. ffT5]) ). Next, consider the case (b) in which the initial value of x^ is small, such as the 
dashed line in Fig. [2] . Since x^ < a/ P^y^b , so x^b will increase until x^b = \j P^>y<f>b = %b ■ 
We have observed above that the y^ evolutionary tracks is not sensitive to the initial value 
of x^ in these two case. But this doesn't mean that xa = %b , since the corresponding N 
might be different, and y^ vary with N . In order to show these two case will converge, we 
suppose that x^, a = x^b + 5 at certain point N (e.g. N = —4). Since x a( p, x^b are very small, 
and in some sense y$ evolutionary track is independent, so the Eq. (fl5l) can be expanded 
to lowest order in 8 : 5' = —35 . The solution of this equation is 5 oc e~ 3N .This means that 
S dacays exponentially, and the evolutionary track of different initial value of x^ given in 
Fig. [2] converge. 

Similar analysis can apply to x v . So we have proved that the evolutionary tracks of 
the Xff,, y^, x^ and y v of different initial value given in Fig.[TJare converge. 

C. The future of the universe 

The critical points correspond to the fixed points where x'^ = 0, y'^ = 0, x' = 0, y' v = 0, 
which have been calculated and given in Table I, and there are self-similar solutions with 

H 3 1 y<i>(l-xl) y^ + x 2 ) \ 

W^-^-^-^TTW--'} <21) 

This corresponds to an expanding universe with a scale factor a(t) given by a oc t p , where 



o / y<t>{~t- x %) _ y v {i+ x %) , \ 

We now study the stability around the critical points given in Table I. Consider small 
perturbations Sx^, 5y^, 6x v , and bx^ about the critical points [x^y^x^y^: x^ c — > 

~t~ o'Xq . y§c * Ufa ~t~ ^y<i>7 x^p C x^p C ~\~ Sx<p 7 y^pc * Vipc ~t~ &y<p' 

Substituting into Eqs. (|T5|) — (|T8|) . lead to the first-order differential equations: 



( 5x '<t> \ 




' Sx^ ^ 




= M 




K 













(23) 



where M. is a matrix that depends upon x$ Cl y^ x^ c and y^. 
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Label 




y<pc 




UipC 


fi)qr>j QfpTl CP 


A. 














all 7 


R 


n 


n 


-.//? 7/ 

V Hipyipc 


2 


51 1 1 '"V 


c. 


±i 











all 7 


D. 


±i 





— a/ QinVmr 


v //3 2 ,+4+ / 3 v 
2 


all 7 


£. 


l 


1 

8a, 








7=1 


F. 


-l 






2 


7 = 1 


G. 




7 








7 < 1(^^ + 4-^) 


H. 




2 








all 7 



TABLE I: The list of the critical points. 



Label 


mi 


m 2 


m 3 


777-4 


Stability 


A. 


-3 


37 


-3 


37 


unstable 


B. 


-3 




— 1(2 + x 2 c ) 


-3(7 + 4c) 


stable 


C. 


6 


37 


-3 


3 7 


unstable 


D. 


6 




— 1(2 + x 2 c ) 


-3(7 + 4c) 


unstable 


E. 





3 
2 


-3 


3 


unstable 


F. 


6 + 6x2 c 


3 T 2 
2 <pc 


— 1(2 + X 2 C ) 


-3(7 + 4c) 


unstable 


G. 


a 


b 


-3 


3 7 


unstable 


H. 


-3 7 + 3x2 c 


-3+^ 


-3 


Q T 2 


unstable 



where: a, b = ^ (^(3 7 - 6) ± 3^16^ 7 2 y/T^j + g(4 - 20 7 + 17 7 2 ) ) 

TABLE II: The eigenvalues and stability of the critical points. 

The general solution for the evolution of linear perturbations can be written as 

5x<f, = Mn exp (miN) + uyi exp {m 2 N) + ui 3 exp (m 3 N) + uu exp {m±N) 
8y<t> = u 2 i exp (m 1 N) + u 22 exp (m 2 N) + u 23 exp (m 3 N) + u 24 exp (m 4 N) 

(24) 

Sxp = u 3 i exp (miN) + u 32 exp (m 2 N) + u 33 exp (m 3 N) + u 34 exp (m 4 N) 
5y<p = U41 exp (miN) + u i2 exp {m 2 N) + u 43 exp {m 3 N) + w 44 exp (m 4 iV) 

where mi, m 2 , m 3 , and m 4 are the eigenvalues of the matrix A4. Thus stability requires 
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FIG. 3: Evolution of the equation of state (cj)and density parameters'^) as a function of 
N for the dark energy model with = /?„ = 1/3, 7 = 1. Initial conditions (at N = —8) 
: a. solid line: = 0.9999999, ^ = 6 x 1CT 11 , x^i = 2.0, y v i = 6.5 x 10~ n ; b. dashed 
line: x^ = 1 x 10~ 12 , y^ = 6 x 10~ n , x^ = lx 10 -12 , y^j = 6.5 x 10~ n ; c. dotted line: = 
0.5, y c / > i = 6x 10~ n , x vi = 1.0, y^i = 6.5 x 10~ n . 



the real part of all eigenvalues being negatives [26j, |29 . 

We obtain the eigenvalues and stability for the fixed points in Table II. The system has 
a fixed point A which is a fluid-dominated solution, a fixed point B which is a phantom 
tachyon-dominated solution, a fixed point C which is tachyon kinetic-dominated solution, 
a fixed point D which is the two-field dominated solution, two fixed points E and F which 
exist only for 7 = 1, a fixed point G in which the energy densities and p 1 decrease with 
the same rate, a fixed H which is tachyon dominated solution. 

In Fig. [3], we plot the evolution of the equation of state in the future. We find that the 
cosmic evolutionary track is towards the stable fixed point B in this model. This can be 
seen by compare Table I with Fig.[2j 

Substituting f3 v = 1/3 into the stable fixed point B ( x^ c = 0, y$ c = 0, x^ c = — y/^/^, 
y vc = (y/P£+4 + P v )/2 ) , we get x^ = 0, V<t>c = 0, X(pc = -0.627285, y vc = 1.180460, which 
are consistent with Fig.[2j And from Eq. (120]) . we know uj = —(1 +x 2 c ) = —1.393487 at the 
fix point B, which is consistent with Fig.0 

Generally speaking, if the initial values of x^, y^, x v and y^ are not the values of the 
unstable critical point given in Table I, they will evolve towards the stable critical point 
B ( If the physical constraints 1 — x| > is not violated ). This point can be seen from 
Eq. ( l24l) . (If one evolutionary track is towards the unstable critical point, when the values 
of x^.y^^x^ and y v differ from x^y^x^y^ by an amount 5, then from Eq. ( 1241) . we 
know the 5 will be larger, instead of becoming smaller. ) 
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D. Discussions 



In a tachyon dark energy model, the tachyon is the only source of the dark energy 



and there are three kind stable critical points [22], which existence depend on the 
value of 7. But in our model, there is only one stable critical point, which existence is 
independent of the value of 7. In a tachyon dark energy model, the value of x^ and 
can be non-zero at the stable critical points. But in our model the value of x^ and y^ 
must be zero, which we have shown in table I. This's not accident. This point can be 
seen as follows. At the critical point, the value of x^ c , y^ c , x vc and y vc are fixed, and 
the value of x^c is non-zero, otherwise the value of y^c is zero (see Eq. (fTTI) ). which 
means p vc = 0, so x vc = is impossible. 

. _ VW _UQ* a , ^ = ». = _ 1 + ^ = _l + 4c >_l (25) 



'c 



2 /I _ r 2 

x 4>c 



V{<p c ) 3M y vc Pipe .2 1 2^-1 

p yc = 7; = = j -z- H u vc = — = -l- Vc =-l- Xipc< -l (26) 

If ?/</>c 7^ 0, then from Eq. ff25l) . we know if 2 is nonincreasing at the fix points, since p$ c 
is nonincreasing. If y 9C 7^ 0, then from Eq. ( |26l) . we know iY 2 is increasing at the fix 
points, since p^ c is increasing. So either y^ c or y vc will be zero . Since p^ is increasing 
and p^ is nonincreasing, so we choose y^ c = 0. And from Eq. (fl5l) . we know = 0. 
(This is because y$ c = 0, and if the physical constraints 1 — x\ c > is not violated .) 

In a tachyon dark energy model with inverse square potential V{4>) = M^<fr~ 2 , in 
order to have significant acceleration at late times(a oc t p ,p = ~ (tjpj ^ 1)j we 



clearly require much larger than the Planck mass 23]. Such a large mass is 
problematic as we expect general relativity itself to break down in such a regime. This 
problem is fortunately alleviated for the inverse power-law potential V = M^~ n 0~™ 
with < n < 2. In our model, since there is another field if, with the equation of 
state uj^ < —1, a significant acceleration at late times can be obtained much easier. 
But the value of still can't be small very much , since if = ^-rA is larger , the 

risk of 1 — <p 2 becoming nonpositive is increasing, which can be see from Eqs. (IT5|) and 

(USD. 



The speed of sound describe the evolution of small perturbations. In a tachyon dark 
energy model the sound speed is 

c 2 = P^L = 1 _ 02 (27) 

P<pX 

where X denotes the partial derivative with respect to X = ^(d^(f)) 2 . Since the value 
of 1 — (p 2 is necessarily nonpositive because of the square root in the Lagrangian density 
Eq. (jSJ), the energy and pressure are real, and inhomogeneous perturbations have a 
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positive sound speed, so the theory is stable. In our model, there are two fields. 
Physically we can use the independent sound speed of each component to describe the 
whole system. However, the present constraints on the sound speed of dark energy are 
so weak that considering to study the two independent sound speed is not justified at 



present. So it can use the effective sound speed, as some authors do [30|, [3jJ. As we 



have shown above, when N is large enough, the fractional energy density of phantom 
tachyon tt v — > 1. So ultimately the effective Lagrangian density is Eq. (EJ), and the 
effective sound speed is 

c 2 = ^X =1+02>1 (2g) 

P^pX 

This means that perturbations of the background scalar field can travel faster than 
light as measured in the preferred frame where the background field is homogeneous. 
But there is no violation of causality. The theory of the k-essence-like scalar fields 
with the Lorentz invariant action is not possible create closed time-like curves in the 
Friedmann universe and hence we cannot send the signal to our own past using the 
superluminal signals build out of the superluminal scalar field perturbations j32j ]. 



III. ANOTHER TACHYON-QUINTOM MODEL INCLUDING THE 
INTERACTION BETWEEN TWO FIELDS 

In order to show some impact of interactions between the two scalars on the evolution of 
the universe, we consider another system, which include a fluid with barotropic equation of 
state p 7 = (7 — l)p 7 , < 7 < 2, and two scalars with interaction between them. Maybe there 
are many more interactions between the two scalars, but for simplicity we only consider the 
below Lagrangian density, to see whether there are some interesting results or not. 

The Lagrangian density of the scalars we choose are: 



C = -V(4>, y?) V 1 + g^^M4> - g^d^d uV (29) 

In this section, we turn our attention to the possibility of the scalars as a source of the 
dark energy. We restrict to spatially homogeneous time dependent solutions for which 
(9j0 = diip = 0. Thus the energy densities and the pressure of the fields are 



P= =, P=-V(<f>,<p)yjl- 2 + v 2 (30) 

yfl - 2 + 2 

Here a dot is derivation with respect to synchronous time. The background equations of 
motion are 



= (31) 
= (32) 
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P 7 = -3n/Hp^ 



(33) 



H 



j?V(<f>,<p) 



<f?v(<f>,<p) 



2M p 



+ 7P 7 



!> 2 + 2 



'1- 2 + 2 \ I I 

together with a constraint equation for the Hubble parameter: 



H 2 



3M p 2 



+ p 7 



1-02 + 2 

The potentials we considered are still inverse square potentials: 

V(<t>, <p) = A x M p 2 0- V -1 + A 2 M 2 0" 2 + A 3 M p V 2 
We define the following dimensionless quantities : 



x<i> = 0, y<j> 



X v = <j>, y v 



Pi 



\/3H' v r ' ^ VSH' 3H 2 M 2 
Now the Eqs. (1351) and ( 1341) can be rewrite as follow: 

+ z = VL DE + z 



_ ^w^Vf + A az/ 2 + A 3 z/ 2 



1 - x\ + X2 



(34) 



(35) 



(36) 



(37) 



(38) 



H' 



3 
2 



1-4 + 4 



7 



(39) 



where Qde measure the dark energy density as a fraction of the critical density, a prime 
denotes a derivative with respect to the logarithm of the scale factor, N = In a. 
Then the evolution Eqs. (l3Tj) and fl32|) can be written to an autonomous system: 



x„ 



-3(1 -xl 



3Ai+3A 2 3/0/y ¥ .+3A 3 y ¥ ,/^ J + X tpJ 



- V3X i - 2y/3 A 3 ^ / 
3Ai+3A 2 ^/y^+3A3y< / ,/y0 



) (i - 4) 



(40) 



y<t> = y4> 



-3Ai^^(7-x|+a:^,)-3A2^(7-3:^,+a:y)-3A3^(7-x|+x^,) 



v^x^ 2 + (41) 
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c: 




FIG. 4: Evolution of the equation of state (oj)and density parameters'^) as a function of 
N for the dark energy model with Ai = 1, 7 = 1. Initial conditions (at N = —8) : a. solid 
line: Xfc = 1000.0, y^ = lxlO -5 , x V i = 1000.0, y v i = lxlO -5 ; b. dashed line: x^ = lxl0~ 12 , y^ = 
1 x 10~ 5 , x^i = 1 x 10~ 12 , y vi = 1 x 10~ 5 ; c. dotted line: = 0.9999999, y^ = 1 X 10" 5 , x vi = 
1 x 10- 12 , y ipi = lx 10~ 5 . 



x' 



-3(1 + x 



-V3A 



-2V%\ 3 yl/ y<t> 



x^ + 



3Ai+3A 2 j/^>/j/ ¥ ,+3A3j/ ¥ ,/i/^ 
- %/3 A 1 j/0 - 2 V3 A 2 /y v 
3Ai+3A2j/^/j/ ¥ ,+3A3j/ ip /j/0 



, 1 X^ 



x , 



(42) 



-3Ai^^( 7 -a:|+^^)-3A 2 ^(7-x|+x^)-3A3^(7--3 : j,+a:^,) 

The equation of state of the dark energy is 

V 



P 



V3x^yl + \ 1Vv (43) 



(44) 



For simplicity, we confine ourselves to the case Ai 7^ 0, A 2 = A 3 = 0, then the Eqs. fj40l — 
fj4~3]) can be reduced to 



x,, 



x 



-3(1 - xl){x^ - -j|(l + x 2 )) - 3x^(x^ + ^|(1 - x\)) 

, -3A l2 / 2 ^( 7 - x 2 + x 2 ) /- 3 

1/0 = 1 v 3x ^ + -7^ 

2 v /l-x 2 +x 2 2 

-3(1 + x\)(x^ + ^|(1 - xj)) + 3x^(x - -J=(l + x 2 ,)) 



(45) 
(46) 
(47) 
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Label x^ 


V4> 




y v 


Existence 


A. 











all 7 


B 4/2 

V 2 


/Iz 
V 2 








all 7 



TABLE III: The list of the critical points. 



Label m\ 


m 2 


m 3 


7714 


Stability 


A. -3 


3 7 
2 


-3 


37 
2 


unstable 


B. i(-3 + 3 v / l + 2 7 2 -4 7 ) i( 


-3 - 3*/l + 27 2 


- 4 7 ) -3 


33 
2 


unstable 



TABLE IV: The eigenvalues and stability of the critical points. 



/ -3An/^( 7 - x\ + x%) R 3 
^ = 1 — - v 3x^7/ + -7^ (48) 

The evolution of the equation of state u and density parameters are shown in Fig. H] with 
Ai = 1, 7 = 1. We choose N = —8 as the initial number of e- folds, so choose the 7 = 1 
in Eqs. (1461) and (j48l is a good approximation. From Fig.H] we can see that this model 
is not sensitive to the initial kinetic energy density of the two fields {x$ = 0, x v = ip). 
From Eq. (130]) and Fig.H], we know that the initial energy density of <fi varied by about four 
orders of magnitude is still consistent with current observational constraints. But the initial 
potential energy density of the two fields require fine-tuning to agree with observations. 

The critical points correspond to the fixed points where x'^ = 0, y'^ = 0, x' = 0, y' v = 0, 
which have been calculated and given in Table III. To study the stability of the critical 
points, we substitute the linear perturbations about the critical points into Eq. (T4"5]) — (148]) 
and keep terms to the first-order in the perturbations. The four perturbation equations give 
rise to four eigenvalues. The stability requires the real part of all eigenvalues be negative (see 
Table IV for the eigenvalues of perturbation equations and the stability of critical points). 

So there are no stable critical points. Since the two kind critical points all have y v = 0, 
which means that z = l(see Eq. (138]) ). The physical constraints 1 — 2 + (p 2 > set limit on 
the equation of state ui (see Eq. (H4"]) ) : uj < 0. Compare this with 7, (7 = 4/3 for radiation, 
7 = 1 for dust matter), we know that the energy densities of the fields p decrease more 
slowly than p 7 , so z = 1 and y^ = are impossible. 



IV. SUMMARY 

In this paper we have studied the tachyon-quintom dark energy models, in which during 
the evolution of the universe the equation of state w changes from w > —1 to w < — 1. 
Firstly, the model we studied is made up of two fields, one is tachyon, the other is phantom 
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tachyon. In order to construct a autonomous system, the potentials we choose are inverse 
square potentials. We find the model is not sensitive to the initial kinetic energy density of 
tachyon and phantom tachyon, and we analyze the reason in detail. The initial energy density 
(at N=8 ) of the tachyon varied by nearly four orders of magnitude is still consistent with 
current observational constraints. The phase-space analysis of the spatially flat FRW models 
shows that there exist a unique stable critical point, and we compare it with tachyon model 
at last. Then we consider another form of two-field model which include the interaction 
between two fields. For the case of Ai 7^ 0, A 2 = 0, A 3 = 0, the phase-space analysis shows 
that there is no stable critical point. In some sense, this work means that multiple kessence- 
like fields can implement the quintom, which extends the possibilities that the quintom is 
realized and is worth further study. 
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